FUNCTION CIRCLE, xcenter, ycenter, radii, density
	num = n_elements(xcenter)
   points = (2 * !PI / (density-1)) * FINDGEN(density)
   x = rebin(xcenter, num, density) + rebin(radii,num,density) * rebin(COS(points), density, num)

   y = rebin(ycenter, num, density) + rebin(radii,num,density) * rebin(SIN(points), density, num)
   RETURN, TRANSPOSE([[[x]],[[y]]])
 END